Collective modes in the electronic polarization of double-layer systems in the 

superconducting state 
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Standard weak coupling methods are used to study collective modes in the superconducting state 
of a double-layer system with intralayer and interlayer interaction, as well as a Josephson-type 
coupling and single particle hopping between the layers by calculating the electronic polarization 
function perpendicular to the layers. New analytical results are derived for the mode frequencies 
corresponding to fluctuations of the relative phase and amplitude of the layer order parameters in the 
case of interlayer pairing and finite hopping t. A new effect is found for finite ^-dependent hopping: 
then the amplitude and phase fluctuations are coupled. Therefore two collective modes may appear 
in the dynamical c-axis conductivity below the threshold energy for breaking Cooper pairs. With 
help of numerical calculations we investigate the temperature dependence of the collective modes 
and show how a plasmon corresponding to charge fluctuations between the layers evolves in the 
normal state. 

PACS numbers: 74.20.Fg, 71.45. Gm, 74. 80. Dm, 74.50.+r, 74.25.Nf 



I. INTRODUCTION 

The problem of collective modes in superconductors-], 
a very old one. Already Bogoliubov and Andersoro 
pointed out that charge oscillations can couple to os- 
cillations of the phase of the superconducting order pa- 
rameter via the pairing interaction. In a neutral system 
this would lead to a sound-like collective mode. In a 
charged system the frequency of this mode is pushed up 
to the plasjna frequency due to the long-range Coulomb 
interactionaa, and at these high frequencies this mode is 
of no importance for the superconducting properties. 

The situation is different for modes which dot-not cou- 
ple to long-range density fluctuations. Leggettcl showed 
that in a two-band superconductor oscillations in the 
occupation difference between the two bands couple to 
phase fluctuations of the order parameters for the two 
bands giving rise to a collective mode with a frequency 
below 2A, the threshold energy for breaking Cooper 
pairs. Also oscillations of the amplitude of the order pa- 
rameter are not perturbed by charge fluctuations. The 
frequency of this mode, however, is at the threshold for 
particle-hole excitations in normal isotropic superconduc- 
tors. Only in special cases overdamping of this mode can 
be avoided!]. Low-frequency collective modes may also 
exist in superconductors with a multi-component order 
parameterD'EI. Such order parameters are frequently dis- 
cussed candidates for some heavy fermion superconduc- 
tors. 

Another possibility to avoid the influence of long-range 
Coulomb forces on the collective modes is realized in 
strongly anisotropic superconductors^], in -particular in a 
periodic system of superconducting layera!3. In the case 
of a finite value of the wave-vector q± perpendicular to 
the layers density fluctuations within the layers do not 
build up long-range Coulomb forces. In the limit q — > 
only the Coulomb interaction between the layers remains. 

Recently the question of collective modes has been 



brought up again in connection with the multi-layer 
structure of the high-T c superconductors. In a double- 
layer system like BSCCO or YBCO one obtains two 
electronic bands for the motion of electrons parallel to 
the layers corresponding to states with symmetric or an- 
tisymmetric wave functions!!!!. In the superconducting 
state both bands acquire gaps. Ih^sugIi systems one can 
discuss different collective modesl!3l!£l corresponding to 
fluctuations in the occupation number of the two bands 
and to charge oscillatiop-|between the two layers. 

In a series of papers! 1 5 H 16 1 we have investigated in de- 
tail the electronic polarization of double-layer systems in 
the normal and superconducting state and have studied 
the influence of charge fluctuations between the layers on 
the renormalization of transverse c-axis phonons. Here 
we have assumed a k-dependent tight-binding coupling 
between the layers and pairing interactions for electrons 
within the same-.and in different layers. In our numeri- 
cal calculations!!!?! we included both the vertex corrections 
due to the BCS-interaction and the Coulomb interaction 
between the layers, thus taking into account the effect of 
possible collective modes. We found a shift of phonon 
frequencies in the superconducting state which is in rea- 
sonable agreement with FIR-experiments. As in this case 
we assumed a fairly large value for the tight-binding cou- 
pling the collective modes are strongly damped by quasi- 
particle excitations. In the present paper we analyze in 
more detail the conditions for undamped collective modes 
in the electronic polarization. The results depend cru- 
cially on the values of the different pairing interactions. 

In principle one has to distinguish three different types 
of interactions: 1. an interaction between electronic den- 
sities within one layer, 2. an interaction between elec- 
tronic densities in different layers. 3. an interaction be- 
tween mixed densities from the two layers. All three in- 
teractions can be mediated by phonons (or other types 
of bosons): the first two by a change of local poten- 
tials induced by lattice displacements, the third one by 
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a change of the hopping energy between the two layers. 
The first two interactions conserve the number of par- 
ticles in the two layers while the third interaction may 
interchange particles between the two layers. In particu- 
lar, this interaction allows a transfer of two particles like 
in a Joscphson coupling. This type of coupling, which 
can also be derived from a second order hopping pro- 
cess has gained much interest recently. P. W. Ander- 
son and coworkers have argued that in strongly corre- 
lated electron systems coherent single-particle hopping 
between the two layers is suppressed, while a coherent 
momentum conserving second-order hopping process is 
possibleEH. In such an interlayer-tunneling model with- 
out single-particle hopping two collective modes involving 
fluctuations of the phase and amplitude difference of the 
layer order-parameters are founded. 

In a two-layer system one has to consider two order 
parameters corresponding to pairing of electrons in one 
layer and in different layers. This leads in general to dif- 
ferent order parameters (and gaps) for the two bands 
which may even have different signs for s-wave pair- 
ing. Depending on the relative sign of the band order- 
parameters one can distinguish between two different 
pairing types, which both have s-wave symmetry: In the 
case of dominant intralayer interaction the order parame- 
ter of the two bands have equal sign, while for .dominant 
interJaiier interaction they have opposite sigrJIaEll. In 
Ref.EjO it is shown that antiferromagnetic interactions 
between the two layers favour a superconducting state 
with interlayer pairing and anisotropic s-wave symmetry. 

In the following discussion we will consider all three 
types of pairing interactions mentioned above and will 
also take into account single-particle hopping between the 
layers. In this paper we consider only s-wave pairing by 
neglecting the k-dependence of the pairing interaction. 
We discuss the collective modes for both pairing types 
and calculate the electronic polarization between the two 
layers, because this quantity enters directly the dynam- 
ical conductivity for electric field vectors in c-direction. 
Our results extend earlier work on collective modes in 
two-layer systems mainly in the following respects: 1. we 
calculate collective modes in the case of interlayer pair- 
ing, 2. we find that a k-dependent tight-binding hopping 
between the layers couples amplitude and phase modes, 
making it possible that both modes appear below the 
threshold frequency in the optical conductivity. 

Our paper is organized as follows: In the following 
section we specify our model and write down the inter- 
actions in a 4 x 4-Nambu matrix notation. Then we cal- 
culate the self-energies and discuss the self-consistency 
equations. In section III we solve the vertex equations 
for the polarization function in the neutral and charged 
system. The collective modes for k-independent hopping 
matrix element are studied in section IV for intralayer 
pairing and interlayer pairing. In section V numerical 
results for the polarization function are presented. The 
results are summarized in section VI. 



II. NAMBU FORMALISM FOR TWO-BAND 
SYSTEMS 

A. Model 

We consider an electronic double-layer system de- 
scribed by the Hamiltonian 

Ho = ^2e k (c\ ka c lklJ + c\ ka c 2ka ) 

ka 

+ t k( C lka C lka+ C lka C 2ka) W 

and the interaction 

— 2 ^11 C lk+qa C ik'-q(j' C ik'<j' C ikcr 

kk'qcrtr' i 

+ VI c L+gcr C ]fe' -qo' C jk' <y' C ika 
+ 3 ( C lk+qa C lk'-qcr' C jk'a' C jk<j 

+ C ik+qa C jk' -qa' C ik' a' C jk<j) (2) 

Here t k describes a tight-binding coupling between the 
two layers i — (1,2), j = 3 — i. The couplings V\\, Vj_, 
J are effective pairing interactions. Vn describes the in- 
teraction between electronic densities within one layer, 
V± between different layers, and J is the coupling be- 
tween mixed densities. The origin of these interactions 
is a combination of an attractive interaction due to the 
exchange of phonons (or other bosons) and a screened 
repulsive Coulomb interaction. Because J includes in- 
trinsic tunneling of Cooper paris from one layer to the 
other, J is called Josephson coupling. Here we neglect 
any momentum dependence of these interactions (except 
of a cut-off introduced later). 

The Hamiltonian Hq can be diagonalized by introduc- 
ing new fermionic operators a aka 

1 , . 1 , , 

a \ka — -J=(C2kcr ~ Clka), O-lka — ~^^ C2k<J + C lka) 

corresponding to states with antisymmetric and symmet- 
ric wave function on the two layers (this symmetry is not 
broken by introducing the interactions). We then obtain 
two bands with quasi-particle energies 

£lfc = £fe — tk, £2k = Cfe + tfc 

In order to treat superconductivity it is usefjul to com- 
bine the Fermi operators to a Nambu spinoiou 

* fe = ( a lfc|' a l-fel' a 2fc|' a 2-fcl) 

With help of these spinors we can express the Hamilto- 
nian of the two-layer system as follows (apart from con- 
stants): 
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(3) 



k,k' ,q 



(4) 



Here V = (V\\ + V x )/2, V = (V\\- V±)/2. D^ 1 are 4 x 4 
matrices in Nambu space which are constructed from two 
sets of Pauli matrices D J = Tjcri. Example: 



D 



13 



a 3 
<J 3 



The terms containing the matrices D 03 , D 33 describe the 
sum and difference of density operators of the two bands, 
while the term with D 13 contains fermionic operators 
from different bands. Accordingly V leads to interband 
transitions while V and J cause intraband transitions 
only. In principle one can derive a further interaction 
mediated by phonons which couples operators character- 
ized by the matrices D 33 and D 03 . Such an interaction is 
omitted here, because it has no influence on the interlayer 
polarization function. 

For the interactions V,V,J a cut-off has to be intro- 
duced either in momentum space or frequency space. We 
introduce the cut-off in momentum space: 



V, V, J 



^ if \e k - \e k < - A*l < l °c 
= otherwise 



(5) 



lu c is the cut-off energy, \x the chemical potential. Note 
that the cut-off has been introduced for the energy e k 
and not for the band-energies e ak . Otherwise we obtain 
inconsistencies in solving the vertex equations. 

In the following we are primarily interested in the cal- 
culation of the interlayer polarization function which is 
the correlation function of the electronic polarization per- 
pendicular to the layers. The latter is described (up to 
a factor ed/2, where d is the layer separation and e the 
electronic charge) by the operator 



P C 2kcr C 2ka 



C lka C lka 



(0) 



ko 



There are three other operators which couple to the in- 
terlayer polarization in the vertex equations. These are: 

^ = — * C 2fcT C 2-fcl — C lfcT C l-fcl _ C 2-fc| C 2 



^ — c 2feT c 2-fei c ifcT c i— fei + c : 



2-fcl C 2fcT C l-fei C lfcT 



Et _ t 

C 2ka C lka C lka C 2ka 



(7) 



In Nambu space a 4 x 4 matrix corresponds to each op- 
erator: 



P jl = ^2^iD ji ^ k 



(8) 



where P = F 13 ,$ = P 12 ,A = P n ,j = P 20 . 

$ and A am. called phase and amplitude operator in 
the literaturet 3 ], because in the case that the phases of 
the layer order-parameters are zero $ and A measure the 
difference of the phases and amplitudes of the layer order- 
parameters, j is related to the interlayer current density. 
This is strictly so only in the case of constant interlayer 
hopping t k — t and vanishing Josephson-coupling J; then 
the equation of motion 2tj — P is fulfilled. In the case 
J ^ the current-density operator becomes more com- 
plicated, it contains also two-particle operators. Never- 
theless j will be called current operator in the following. 



B. Self-consistency equations 

The Green's functions can be combined to a 4 x 4- 
Nambu matrix. We assume the Green's function matrix 
to be diagonal in the band indices: 



G(k,z) 



Gi(fc,z) 
G 2 (k,z) 



(9) 



where G a (k,z) are 2x2 matrices for each band. This 
assumption which neglects pairing in different bands is 
reasonable since the two bands (in the presence of aj£,- 
nite tight-binding coupling) have different symmetriesEil, 
and this symmetry is preserved by the interaction H$. 
One consequence of this assumption is, that the phases 
of the order parameters < Ci k ^Ci^ k ^> for pairing in one 
layer are equal for the two layers (equal amplitude and 
phase). Intraband pairing seems to be the ground state 
for all negative values of J, while it may become un- 
stable for large positive J 3> \t k \. This is supported 
by the calculation of the groundstate of a model with 
two electrons on two sites with all the interactions con- 
tained in the Hamiltonian Hq + H$ ■ Defining operators 
a aa = ^(c2 CT ±cife .) for the single-particle eigenstates of 
Hq one finds that the groundstate contains only diagonal 
terms a} a ^a} a ^ |0 > as long as the Josephson coupling is at- 
tractive J < but contains non-diagonal terms for large 
positive values J> \t k \. In our numerical calculations of 
the polarization function we will show that an instability 
occurs for large positive Josephson coupling J. 
The bare Green's function is given by: 



G 0a( k > z ) = Z(J o - (e«,fc - A*) ^ 



(10) 



In the presence of the pairing interaction we have self- 
energy corrections: 



ka 



G- 1 (k,z) = Go 1 {k,z)-E(z) 



(11) 
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where the self-energy contains contributions from 

intraband and interband interactions. 

E = ~ VD 03 G(k, iuJ m )D 03 + VD 13 G(k, iuj m )D 13 



JD G(k, ito m )D 



(12) 



The self-energy is diagonal in the band indices 

^(» e") < 13 > 

and the 2x2 components S a are given by 

51 = -- V" Wa 3 Gi(k,iu> m )a 3 + Va 3 G 2 (k,iu> m )<j 3 

P i 

k.Lj m 

5 2 = V Wa 3 G 2 (k,i(J m )a 3 + V<7 3 G 1 (k,ioj m )a 3 (14) 

ku m 

with the intraband interaction W = V + J. While the 
interactions V, J enter in the same way into the self- 
energy, this is not so in the vertex equations. Therefore 
J cannot simply be incorporated into the interaction V 
by a proper redefinition. 

Writing E Q = — A a ai — x a a 3 we have 



G a 1 (k, z) = zero - £, ak + A a ai 



(15) 



with £ Qfc = e ak +x a - fi. 

A a is the energy gap of band a. x a is an energy shift 
due to the pairing interaction which is also present in 
the normal state. In a one-layer system be in- 

corporated into the chemical potential /x, but this is not 
possible in the two-layer system with interlayer hopping. 

Performing the frequency summation we get the fol- 
lowing self-consistency equations for the gaps and energy 
shifts: 



Ai = -WAih - VA 2 I 2 
A 2 = -WA 2 I 2 - VAxh 

x x = -WKi - VK 2 
x 2 = -WK 2 - VKi 



with the integrals 

El j ]_ 0E a k 
OR. 



I n = 



■ tanh : 
2E n k 2 



(16) 



(17) 



Cafe , n (3E a k 

■ tanh 



2E n 



and E ak = y/£ k + A| 

Without a cut-off the integrals Ii would diverge. Ap- 
proximately the integrals have the values: 



I a ~ N a \n(2uj c /Max(A a ,T)) 



(18) 



where N a is the density of states of band a per spin 
at the Fermi surface, and iv r is the cutoff. 



The self-consistency equations (|l7|) describing the en- 
ergy shifts X1.2 of the two bands can be solved ap- 
proximately by observing that the integrands are nearly 
step functions. Assuming constant density of states N±, 
N 2 (this is given by a quadratic 2D-dispersion) and k- 
independent interlayer hopping t one finds 



xi 
x 2 



-WN 2 {t~x 2 ) 



VN 2 {t-x 2 ) 
+ VN 1 (t + x 1 ) 



(19) 



In the limit of W, V — * and cut-off lu c — > 00 these 
shifts vanish, but the ratio (x 2 — X \)/(W — V) stays finite 
(N - (JVi+JV 2 )/2): 



x 2 - Xl 
W -V 



-2N t 



(20) 



The solution of the gap equations_for a two-band sys- 
tem is discussed already by LeggettEl; the qualitative re- 
sults can be summarized as follows: There exists no so- 
lution, if W > and W 2 > V 2 > 0; then the intralayer 
TTj = V]| + J and the interlayer coupling W± = V± + J are 
repulsive. There are two nontrivial solutions, if W < 
and W 2 > V 2 > 0; in this case W11 and W± are attractive. 
The solution corresponding to the ground state depends 
on the sign of the interband interaction V. If V < 0, 
the state where the gaps have the same sign (A1A2 > 0) 
is stable; if V > 0, then the gaps have opposite signs 
(A1A2 < 0). The reason is a term in the free energy pro- 
portional to l^cos(</>i — <p 2 ), where 4>i,4> 2 are the phases 
of the order parameters Ai, A2. 

If the density of states for the two bands are equal 
(this happens for a constant density of states and for 
independent interlayer hopping t) the two integrals I t are 
equal for equal values of the gaps (I\ = I 2 for A^ = A|) 
then for all solutions of the gap equations one always finds 
|Ai| = |A 2 | and x\ = — x 2 irrespective of the values of 
the coupling constants, i.e. one either has pure intralayer 
pairing for W\\ < with (Ai + A 2 )/2 = A|| ^ 0, (A 2 - 
Ai)/2 = Aj^ = or pure interlayer pairing for W± < 
with Aj^ 7^ and A|| =0. In the more general case the 
density of states of the two bands are different at the 
Fermi surface, and the two integrals have different values 
also in the case of equal or vanishing gaps, and the two 
order parameters are coupled. 

For the numerical solutions of the self-consistency 
equations and vertex equations with k-dependent hop- 
ping we use the following simple model by choosing two 
bands with quadratic dispersion but different effective 
masses (and different density of states) : 



eifc 

£2fc 



Eft 



tk 
t k 



k 2 /(2 mi ) - t 
k 2 /(2m 2 )+t 



(21) 



In the numerical calculations we use the following dis- 
persion parameters to get k-dependent hopping: 

2 mi /h 2 = leV^A- 2 , 2m 2 /h 2 = 1.2 eV" 1 !" 2 , 

t = 0, 11 = 0.3 eV, uj c = 0.25 eV (22) 
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For the intraband and interband coupling constants W = 
V + J and V, entering the self-consistency equations we 
choose 

N W = -0.139, N V = ±0.185 (23) 
We have chosen these values for the coupling constants 




|k| (1/A) 

FIG. 1. Dispersions in the normal state e%k, £2fc and 
quasi-particle dispersions in the superconducting state (T=0) 
in the fc-region \ek\ < ijJ d . The parameters (22,23) are used. 



in order to get T c -values which are appropriate for high- 
T c superconductors and ensure also the appearance of 
the collective modes we want to study. 

Fig. [l] shows the dispersion for the two bands in the 
normal state and the quasi-particle dispersions in the su- 
perconducting state at zero temperature in the k-region 
|efe — fj,\ < lo c . Figures |2|(a) , ||(b) show the temperature 
dependence of the gaps for negative V (a) and positive V 
(b). Because the effective masses mi and 1112 are nearly 
the same (P4) the values of the two gaps Ai and A2 are 
nearly the same. For repulsive interband coupling V the 
stable gaps have different signs. In both cases intralayer 
and interlayer pairing are mixed. But for negative V in- 
tralayer pairing dominates, for positive V interlayer pair- 
ing dominates. T c is 155 K. We have also calculated the 
level shifts Xi. They are practically temperature inde- 
pendent. 



III. VERTEX EQUATIONS 

Now we will set up the integral equations for the ver- 
tex functions corresponding to the interlayer polarization 
and related operators defined in (^) for the neutral and 
charged system. We are interested in optical response 
functions for electric field vectors perpendicular to the 
layers. Here we can assume that the external wave vec- 
tor is zero. 



25 



20 



> 
CD 

E 15 



10 



20 



> 

CD 

E 



-20 





(a) 


1 1 1 1 1 1 — 


(b) 















16 



4 8 12 

Temperature (meV) 

FIG. 2. Gaps as function of temperature for negative 
(a) and positive interband coupling V (b). The parameters 
(22,23) are used. 



A. Neutral system 

We consider response functions where the bare vertex 
is one of the matrices D^ 1 , corresponding to the operators 
(^J). In the standard ladder approximation the vertex 
equation for the renormalization of a vertex with matrix 
D^ 1 reads 



T jl (w s ) = D> 1 
-\y VD m G{k,iu n + iu s )Tl l {iu 8 )G{k,iun)D 



0:S 



13 



+ VD 13 G(k,iuj n + iuj s )T jl (iuj s )G{k,iLu n )D 



+ JD 66 G{k,iu n + iuj s )T : > l {iuj s )G{k,iuj n )D 66 (24) 



The vertex function depends only on the external fre- 
quency u> s , not on the momentum k (apart from the cut- 
off) and the internal frequency of the Green's function, 
because the interactions are assumed to be constant. The 
label jl indicates to which bare vertex the vertex function 
belongs. In the case of the polarization we have to con- 
sider the matrix D 13 . Of course, due to the interactions 
the renormalized vertex T^ 1 contains also contributions 
from other matrices. 
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The interactions V, J induce only intraband transi- 
tions, the interaction V only interband transitions. These 
interactions do not change the off-diagonal character of 
the bare vertices, therefore we can make the following 
ansatz for the renormalized vertex of the neutral system 



7 
7 



with 2x2 matrices j(iuj s ), ^{iuj s ). Then the system of 
vertex equations can be written as two coupled equations 
for 2x2 matrices: 



7 + W'B + VB = I 
7 + W'B + VB = I 



(26) 



where W = V-J = (V\\+V±)/2-J. The 2 x 2 matrices I 
and I are the upper right and lower left part of the matrix 
D jl . In the case of the polarization D 3 we have I = 03, 
/ = —03. The quantites B arc defined as 



B(iui s 



-iw s )7(iw s )G 2 (fc, iuj n )a 3 



B(iuj s ) = - ^2^2 CT 3G 2 (fc, iuJ n + iu s )i(iuj s )Gi(k, iuj n )a 3 

'■' k 

(27) 



depending linearly on ^(iui s ). Note, that in contrast to 
the self-consistency equations ( PI |l7|) (W = V + J) the 
interaction J enters here with an other sign (W = V—J). 

For the solution of this system of equations it is con- 
venient to decompose 7,7 and also B,B into Pauli ma- 
trices: 



7 = 7oO"o + 71 °i + 72*0-2 + 730-3 

-B = _B O"0 + Bx(J\ + 52*0-2 + -B3CT3 



(28) 



If we write the coefficient functions 7.; as vector 7 = 
(70 j lii 727 73) w e can express the linear dependence of 
B on these coefficients in matrix form: 



B = Kj, B = K^i 



(29) 



The 4x4 matrices K, K are integrals over products of 
Green's functions. For notational convenience it is useful 
to introduce an extra factor i for the ^-component in 
(p8|). The structure of K is given in the appendix (A3) 
and the functions are listed in (A5). With help of these 
matrices K, K the vertex equations can be written as 



(1 + W'K)2+VK2 = L 
(1 + W'K)i + VKj = I 



(30) 



The matrix K differs from the matrix K only by the 
sign of some elements in the first row and column. More 
precisely: K = gKg with g being the diagonal matrix 



g = diag(— 1, 1, 1, 1). I is a vector given by (0,0,0,1), 
(0,0, -i,0), (0,1,0,0), (-*, 0,0,0) for P 13 , P 12 ,P n , P 20 , 
respectively. I_ fulfills: j = gl. The second equation in 
(p0|) can be reduced to the first equation by setting 



1 = 91 



(25) and it is sufficient to solve the first equation 



J2 + i W ' + V 9ii)Kij) ij = h 



(31) 



(32) 



which is a set of 4 linear equations. 

A closer inspection of the integrals contained in the 
matrix K shows that two of the integrals, K\\ and i^22 
are badly convergent. They depend logarithmically on 
the cut-off in momentum space in a similar way as the 
integrals in the self-consistency equations (|l^). These 
badly c onvergent integrals can be eliminated by using re- 
lations ( B16 . BI7 ) which can be derived from Ward iden- 
tities. Then the cut-off is only needed in the calculation 
of the gaps A< . 

The polarization function of the neutral system is given 

by 

< P;P> iWs = 

\ T r{D 1Z G(k, iuo n + iu s )T 13 {tu s )G(k, iuo n )} (33) 

' Wn k 

The polarization function can be expressed by the ma- 
trix B (|27]) or by the matrix K using the commutation 
relations for the Pauli matrices and the development into 
Pauli matrices 



< P, P >=< P 13 , P 13 >= 4K 3j j 



r? 



(34) 



B. Charged system 

Up to now we have only considered neutral supercon- 
ductors. The long-range Coulomb interaction is knowi. 
have important consequences for the collective modesu 
Here we are interested in the optical properties for 
field vectors perpendicular to the layers in the long- 
wavelength limit q — ► 0. Then we have to consider only 
the Coulomb interaction arising from charge fluctuations 
between the layers. These Coulomb forces stay finite in 
the long-wavelength limit but nevertheless are sizable. 
They can be incorporated in the vertex equations in a 
RPA-type manner. 

Let T ll (iuj s ) be the vertex for the charged system. The 
vertex equation for T l1 is obtained from the vertex of the 
neutral system (p4|), if T 3 ' is replaced by T tl (iuj s ), and on 
the r.h.s. of the vertex equation (^4j) the following term 
is included: 



D 13 v^2Y^ Tr {D 13 G{k,iuj n + iuj s )T jl (ioj s )G(k,iuj n )} 

(35) 



G 



with v = lini| gH0 e 2 7r(l - e~ d W)/(\q\e Q ) = e 2 nd/e the 
Coulomb interaction between the layers. For a periodic 
double-layer system v has to be replaced by ire 2 d(c — 
d)/ceo, where d is the distance between the two layers in 
one unit cell and c the distance between the unit cells. 
The vertex equation for T J can also be expressed by 
correlation functions and the vertex of the neutral system 
r ji , as shown in Fig. ^[ This equation reads: 



T jl = T jl + T i3 5 < pl3 pjl > 



(36) 



where the correlation function <C P 13 , P^ 1 S> of the neu- 
tral system to which the interband Coulomb-interaction 
couples is given by: 



u, k 

Tr{D 13 G{k,iu n + iLu s )T l1 (iu s )G{k,iu n )} (37) 



In the following we need the special form of the matrix 
K given in the Appendix: 



K = 



( H++ 


-Y 


-z 


X 


\ 


Y 


-H-+ 


X 


z 


Z 


X 


-H— 


Y 




K X 


-Z 


-Y 


H+- 


) 



(41) 



Furthermore we need the following combinations of cou- 
pling constants, which are listed here for completeness: 



W = V + J, W\\=V\\ + J, W ± ^V ± -J 
W' = V-J, Wn =V\\ — J, W' X =V±-J 



(42) 



with V = (Vy + V±)/2, V = (V]| + Vj_)/2. Of special 
importance is the limit, where all pairing interactions are 
small 

WV , VN Q , J7V «1 (43) 

This will be called the weak coupling limit in the follow- 
ing. 





+ 



FIG. 3. Diagrammatic representation of the vertex T ji of 
the charged system with help of the vertex r ji (full circle) 
of the neutral system. The interband Coulomb-interaction 
v couples to the correlation function <C P 1S , P S> of the 
neutral system. 



From this the vertex function T 13 for the electronic 
polarization is easily calculated: 



-13 



r 1 7(l -v < P L \p Li >) 



(38) 



and we get for the polarization function for the charged 
system: 



« P, p »£, =« p, p » 4Ws / (l - v « P, P ) 



(39) 



IV. ANALYTICAL DISCUSSION OF THE 
COLLECTIVE MODES 



In general the 4 vertex equations ([32 



(40) 



have to be solved numerically. But in some limits the 
equations can be simplified so that analytical results for 
the collective modes can be derived: 1. in the normal 
state and 2. in the limit of momentum independent tight- 
binding coupling tfc between the layers. 



A. Normal state 

In the normal state the functions Z, Z, Y, Y van- 
ish, because they are proportional to the gaps JA5|). 
Therefore the set of equations (|40|) decouples. With 
I = (0,0,0,1) we obtain for the polarization vertex 

'V — <yt3 _ *y-P; 



7o = -W' ± X/N, 7i = 72 = 0, 
73 = (1 + W' L H ++ )/N 



(44) 



with N = (l + W' L H ++ )(l + W( l H + -)- W^W'^X 2 . The 
polarization function ( p4| ) is then given by 

« P, P >^ s = 4 ((1 + W L H ++ )H+- - W ± (X) 2 ^j /N 

(45) 

In the case of momentum independent hopping (t^ — t) 
this expression can be simplified considerably by using 
the relations between the different matrix elements de- 



rived from the Ward identities ( B15 ) and the approxi- 
mation (#2 — Xi)/W± — —2tNo. We obtain 



< P; P »i 



16N t 2 (l + 2N J)/M 



(46) 



with 



M = (iuj s ) 2 - 4t 2 [l - 2(V - 2J)N - AJ(V - J)N$] 

The polarization function has a simple pole which is 
shifted with respect to the simple particle-hole excita- 
tion energy 2t due to the pairing interactions. In the 
weak coupling limit (^) we recover as result the bare 
polarization function, which is given by . In the 
charged system there exists then a plasmon describing 
interlayer charge fluctuations at uo 2 = (2i) 2 (l + vN ). 
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B. fc-independent hopping matrix-element t 

In the case of momentum independent tight-binding 
coupling (£fc = t) the solution for the self-consistency 
equations (17, 161) has the form 



Ai = ±A2 and x\ 



-X2 



(47) 



independent of the details of the band dispersion e^. 
Then there exists a symmetry relation between the band 
dispersions £ife = £fe — t - x i an d &k = £fc + t — x 2 
= Sk — /J,) and the quasi-particle dispersions: 



Ei(-tk) = Mb) 



(48) 



Due to this symmetry some matrix functions Kij vanish 



and therefore (40) is reduced to a system of three coupled 
equations. We have to distinguish the two cases Ai = 
±A 2 . 



1. Intralayer pairing Ai = A2 

In the case of dominant intralayer interaction W» the 
band gaps are equal (Ai = A2 = Ay) and interlayer pair- 
ing vanishes (Aj_ = (A2 — Ai)/2 = 0). Three functions 
are zero because of the relation (|4S 



Z(iuj s ) — X(ilu s ) = Y(iuj s ) = 



(49) 



Therefore the second equation of the vertex equations 
( pp| ) can be solved easily 

lilt, 7f = 0> 7i A (^) = 1/ (1 - W l{ H-+(iu s )) , 

7 ; A = 0, I = 0, 2, 3 (50) 

This means that oscillations of the amplitude A do not 
couple to oscillations of the current j, polarization P, 
and phase <&. But the density, phase, and current oscilla- 
tions couple to each other. The Coulomb interaction does 
not affect the amplitude fluctuatio:p&.|-4rhis is in agree- 
ment with the results found in RefsJl3'EJ, where a double 
layer system with t = is investigated. In the case of 
k-dependent hopping all these quantities couple, which 
will be demonstrated below by our numerical results. 

The frequencies of the collective modes are determined 
by the zeros of the denominator of the vertex functions 
7 3 '. For the current, phase and polarization oscillations 
the frequency of the collective mode is given by the deter- 
minant of the cor resp onding 3x3 matrix of the vertex- 
equation system (|40|); the collective mode of the ampli- 
tude oscillations is determined by (50). 



The functions K(u)+i5) obtain imaginary parts at tem- 
perature T = 0, when the frequency of the external field 
is large enough to break Cooper pairs (ui > min(Eik + 

E 2k ) = J{2t) 2 + (2A,,) 2 with t = t+{xi- x 2 )/2). We 



are interested in undamped collective modes below the 
particle-hole threshold. In order to derive a simple ap- 
proximation for the collective mode frequencies, we shall 



assume small frequencies to <§; y (2i ) 2 + (2Ay) 2 , a small 
hopping matrix element (t/Au) 2 <t^. 1 and zero tempera- 



ture. 



The functions K (A5) can be evaluated with help of 
an expansion in i/An and w/2Am: 



H+-(uj) « -N + uj 2 C, H++ « (2t) 2 C 
Y(co) ^cu(N a + ((2t) 2 - u 2 )C) /(2A,|), 

Z(u) w — 27Vo/(2A|i), X w -uo2tC 



(51) 



where C = -iV /(8A|). Here -(x 2 - x x )/{W±_2i) has 
been approximated by the density of states at the Fermi 
level N (see ©). 

With these expressions we expand the determinant up 
to second order in t/A» or u>/(2A») 



Det 



N W(, 
(2A0 5 



(l-a; 2 W]C) (oj 2 ~u 2 P ) 



(52) 



Here up is the frequency of the collective phase mode in 
the neutral system 



u? P = [loI + {2tf (1 - W'ANo - u 2 C))] j 



1 - WIiNq 



u 2 W'C 



(53) 



with w„ 



(2A|, 



-2J 



No V 2 - J 2 

This reduces to a simple result in the weak coupling limit 



(2tf 



(54) 



Because we have assumed small frequencies the formula 
is only valid for small wo, i.e. \J/Vn \ has to be small. The 
expression for the frequency u>o has been derived already 
by Leggetra. Here we see how the mode-frequency is 
modified by the interlayer hopping t. It is justified to call 
this mode phase mode, because it appears as a resonance 
in the polarization function due to its coupling to the 
phase fluctuation. 

Evaluating the polarization function (J|J) of the neutral 
system for small frequencies and small hopping matrix 
element we obtain 



< P, P »u 



AN, 



u 2 P {\ - wgWjfCT) - u 2 u 2 C (1 - W(\N ) 
(cj 2 - Lo 2 p ) (1 - ^ 2 W/]C)(1 - Wp ) 

(55) 



Inserting (|5^) into the polarization function with 
Coulomb interaction (^) we get the phase-mode fre- 
quency in the charged system: 
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(2t 



t s i2 



Tyi(7V -^o 2 C))] 
l + JV (4i;- Wjf) 



1 + lo%C(4v-W(,) 



(56) 



In the weak coupling limit 
have 



|) and for small |w 2 C| we 



u%= (loI + {2tf){\ + m v) 



(57) 



The interband Coulomb interaction v shifts the mode up 
to higher frequencies, v is proportional to the distance 
d between the layers. In order to get a phase mode be- 
low the particle hole threshold, it is necessary that the 
distance of the layers and the Josephson coupling con- 
stant J is small. If the Josephson coupling J is zero, the 
resonance appears at 2t as in the normal state but still 
below the particle-hole threshold. The mode corresponds 
to the collective tunneling of Cooper pairs between the 
layers without pair breaking. 

The denominator of 7^ ( [SOl) d etermines the amplitude 
mod e. W ith the relations ( |B18[ ) and using Ward identi- 
ties ( B15| ) we get: 



7i» 



= 1 - 



W \\ 



(2t; 



(2An 



1 Wi 



,2A||F(w) 



(58) 



For vanishing Josephson coupling (W^j/Vt^i = 1) the am- 
plitude mode lies just at the particle-hole threshold: 

u\ = (2t) 2 + (2A||) 2 

For finite Josephson coupling an approximate analytical 
result can be obtained with (|5l|): 



(2£) 2 + 2K 2 + (2A|,) 2 ) 



(59) 



Because we assumed small frequencies, this result is only 



valid, 



is negative, which is obtained for positive 



J. For the general solution we refer to the numerical 
calculation. 

The formula for the phase (|^|) and amplitude mode 
( |59| ) show how the mode frequencies are shifted by the 
hopping t and the results of Refsfl'Ej are modified. Note, 
however, that the particle-hole threshold is shifted too. 



2. Interlay er pairing Ai 



In the case of dominant interlayer interaction, when 
A ± = (A a - AO/2 ^ 0, A,, = (A 2 + A x )/2 = 0, the 
following three functions are zero because of (|48|): 



Z(iuj s ) = X(ilu s ) = Y(iuj s ) = 



(60) 



The third column and line of the vertex equation matrix 
( pp| ) is zero with the exception of the diagonal element. 



Now the oscillations of $ are decoupled from the the os- 
cillations of the other quantities, while j, A and P couple 
to each other. 

Again we study the limit of small frequencies lo <C 



(2t) 2 + Aj_ and small hopping |t/A_L | 2 < 1. The func- 
tions can be approximated as: 



Y 



« -n + uj 2 c, h+- pb (2t yc, 

u(N + ((2tf -lo 2 )C)I (2A X ), 

Z w -2&V /(2Ax), X w -u)2tC, 
H-+ » l/W ± - (lu 2 - (2t) 2 )N Q /(2A ± ) 



(61) 



where C = — JVo/8Aj_ and t = t + (x\ — x%)/2 as before. 
From the vanishing of the determinant of the three cou- 
pled equations we now obtain for the frequency of the 
collective mode determining polarization and amplitude 
oscillations: 



+ (2t) 2 (l-W^No-ulC)) 



1 - W'j_N Q 
oj 2 W'^C 

(62) 



1 



with 



(2A ± ) 2 2(V - J) 



u iV V 2 - (V - J) 2 

In the week coupling limit ( f43|) the mode in the neutral 
system becomes: 



= (2t) 2 +u 2 



(63) 



Expanding the nominator and denominator of the po- 
larization function ( pi|) up to second order in i/ A± and 
cij/2A_l we obtain 



«P,P» tJ =4(7V -^ 2 C) 



(2t) 2 (l~ W^Np) 

(UJ 2 - L0 2 p ) (1 - UJ 2 W'^C) 



(64) 



In this case the polarization function is zero, if the single- 
particle hopping t vanishes. This is in contrast to the 
case of intralayer pairing, where the pair tunneling due 
to J produces density fluctuations between the layers 
and therefore spectral weight in the polarization function 
even in the absence of single-particle hopping. P. W. An- 
derson argues that due to correlations the single-particle 
hopping is suppressed. If this is right, then in the case 
of interlayer pairing no excitation of charge fluctuations 
between the layers should be seen in the optical c-axes 
conductivity cr(u>), which is related to the polarization 
function by: a{uo + iS) = ilu << P, P »u+i8 {ed/2) 2 . 
The collective mode in the charged system is given by: 



f (2t) (l + (4w- 

1 - W' ± N 
\-uj 2 W' ± C 



W(MN - u 2 C) 



(65) 
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Only the term containing the single-particle hopping is 
influenced by the Coulomb interaction. This is consistent 
with the picture of interlayer pairing. The hopping of 
interlayer Cooper-pairs produces no charge fluctuation 
between the layers. When an interlayer pair is tunneling, 
one electron jumps to the upper layer, the other to the 
lower layer and the charge density on the layers does not 
change. 

In a similar way as above we obtain for the frequency 
of the phase mode, which is now decoupled from the po- 
larization: 

c4 = (2£) 2 + 2(cj 2 + (2A ± ) 2 ) (66) 

Because we assumed small frequencies, this result is only 
valid, if ujq is negative, which is obtained for negative 
V — J. For the general solution we refer to the numerical 
calculation. 

In comparison to the case of intralayer pairing the 
modes have changed their role. The amplitude mode 
is now the low lying excitation and is coupled to the 
electronic polarization. As in the ground state of inter- 
layer pairing the layer order-parameters < C\-\C\^ > and 
< C2]C2[ > are zero it is no longer justified to call os- 
cillations of <£> phase oscillations. We kept this name for 
convenience only. 

V. NUMERICAL RESULTS 

In the general case of k-dependent hopping when in- 
tralayer and interlayer pairing are mixed it is not pos- 
sible to derive simple expressions for the dispersions of 
the collective modes. We made numerical calculations for 
the polarization functions with the dispersion parameters 
( p2| ) and coupling constants 

N W = -0.139, N V = ±0.185, N v = 1.3 (67) 

for constant W = V + J and variable Josephson coupling 
J for dominating intralayer (V < 0) and interlayer pair- 
ing (V > 0), in order to study the influence of the Joseph- 
son coupling J and the temperature dependence of the 
collective modes. We present results both for the neutral 
system and taking into account the Coulomb interaction 
v between the layers. To resolve the 5-peaks of the col- 
lective modes at T = below the particle-hole threshold 
we have introduced a small imaginary part (10 -3 meV) 
to the frequency. 

A. Dominant intralayer pairing 

First we discuss the results for the case of dominant 
intralayer pairing (V < 0,AiA 2 > 0). In Fig. | the 
spectrum of the collective modes (sharp spikes) and the 
form of the particle-hole spectrum at T — is shown. 
The imaginary parts of the polarization function in the 



neutral |||(a) and the charged system |](b) are calculated 
for different Josephson couplings J and constant intra- 
band W = V + J and interband coupling V . For com- 
parison results of the approximation formulas are plot- 
ted, too. The rhombs and crosses in (a) refer to the 
formulas for the phase ( |53| ) and amplitude mode (|59|), re- 
spectively, calculated with a constant averaged hopping 
t = — 30.5meV. The dashed line represents the phase 
mode frequency (|54|) in the weak coupling limit. 

For negative Josephson couplings J a collective mode, 
the phase mode, lies just below the particle-hole thresh- 
old. With increasing positive J the phase-mode peak 
moves to lower frequencies. 

For positive J another mode, the amplitude mode, be- 
comes visible near the particle-hole threshold. For in- 
creasing J the peak moves to smaller frequencies. Thus, 
two superconducting collective modes can occur in the 
polarization function perpendicular to the layers. This 
effect is caused by the coupling of the phase and ampli- 
tude fluctuations. 

The two approximation formulas for the phase mode, 
represented by the rhombs and the dashed line, give qual- 
itatively the right position of the collective-mode frequen- 
cies, as long as the peaks are below the particle-hole 
threshold. The same holds for the approximated values 
for the amplitude mode (crosses). 

For No J > 0.06 the peak of the phase mode has passed 
zero, and its frequency has become imaginary. This be- 
haviour is also obtained by using the approximation for- 
mula (|54l). We believe, that this instability indicates 
a phase transition from intraband to interband pairing. 
This is supported by a study of a two-site model with two 
electrons (see section self-consistency equation), where a 
strong positive coupling J leads to a ground state con- 
sisting of interband Cooper-pairs (a^a^). 

In Fig. ^(b) the influence of the Coulomb interaction 
between the layers is shown. A plasmon exists within 
the particle-hole continuum. This mode is caused by 
charge fluctuations between the layers. The damping of 
this mode increases with increasing |«7|. The particle- 
hole threshold peak is strongly suppressed, if there ex- 
ists a collective mode with large spectral weight above 
the threshold. Again the approximation formula for the 
phase mode in the weak coupling limit, represented by 
the dashed line in Fig. ^(b), gives the right position of 
the plasmon peak for small \J\. The plasmon mode and 
phase mode coincide for small |J|. 

But for positive J a peak below the particle-hole 
threshold appears (N J = 0.02, 0.04). This peak can 
be attributed to phase fluctuations by comparing the 
peak positions with the analytical results obtained for 
constant t (rhombs). The amplitude peak lies just at 
the particle-hole threshold at N$J = 0.04 and can be 
seen at N$J = 0.06. Those peaks which we have at- 
tributed to the phase and amplitude modes have large 
spectral weight in the correlations functions <C $, $ ^ 
and <C A,A^S>, respectively. 
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Now we discuss the temperature dependence of the po- 
larization functions for large positive (0.05) and negative 
(-0.06) Josephson coupling. 

Fig. |^ shows the imaginary part of the bare polariza- 
tion function (a), the polarization function in the neutral 
system (b), and the real part of the optical conductivity 
in the charged system (c). The latter is related to 
the polarization function by cr[u> + iS) = iuj(ed/2) 2 <C 
P, P ^ n picture ||(a) only particle-hole excitations 

with energy \E\k + ^2fc| are visible at T — 0. For finite 
temperatures also excitations with \Eik — -£?2fc| from ex- 
cited levels are possible at lower frequencies. The broad 
peak at T c is due to interband transitions with different 
tfc-values. The Fig. |](b) shows the complicated transi- 
tion from two distinct collective modes at T = (phase 
and amplitude mode) to a collective density fluctuation in 
the normal state. The Coulomb interaction between the 
layers shifts the collective modes associated with density 
fluctuations to higher frequencies (Fig. 0(c)), whereas 
the amplitude mode near the particle-hole threshold is 
only slightly shifted. 

In the case of negative J, which is the more realistic 
one, if this coupling is produced by the exchange of a bo- 
son, the Coulomb interaction shifts the collective modes 
into the region of particle-hole excitations. Fig. || shows 
the temperature dependence of the remaining density- 
fluctuation spectrum. 

B. Dominant interlayer pairing 

In the case of dominant interlayer pairing (V > 
0,AiA2 < 0) it is possible to produce the same peak 
structures in the polarization function as in the case of 
dominant intralayer pairing. Fig. ^ shows the imaginary 
parts of the polarization function in the neutral (a) and 
charged (b) system for selected J-values. The structures 
are nearly the same as in Fig. [|. For J greater equal 
V collective modes below the particle-hole threshold can 
occur in the polarization function of the neutral system 
and for relative large positive J also in the polarization 
function for the charged system. The damping behavior 
of the plasmon is quite similar as in the case of dominant 
intralayer pairing. For negative and small positive J the 
plasmon peak is strongly damped. 

An example for the temperature dependence of the 
imaginary part of the polarization function in the neutral 
(full lines) and charged (dashed lines) system shows Fig. 
B. The plasmon (dashed line) has the same behaviour as 
in Fig. g. 

VI. CONCLUSIONS 

In this paper we have studied collective modes in su- 
perconducting double-layer systems which are connected 
with charge fluctuations between the layers. In our model 



we assumed a single-particle tight-binding hopping tk be- 
tween the layers and considered three different types of 
pairing interactions: an intralayer interaction Vm, an in- 
terlayer interaction V± and a Josephson-type coupling J 
which allows a transfer of two particles between the lay- 
ers. With these interactions we set up a system of vertex 
equations within a conserving approximation. The col- 
lective modes then appear as resonances in the vertex 
functions. So we were able to study the interplay of the 
three types of interlayer interactions tk, V±, and J, which 
has not been investigated so far together. 

Due to the interaction V± Cooper pairs with electrons 
in different layers are formed. Depending on the relative 
size of the different interactions one finds dominating in- 
terlayer or intralayer pairing. In the case of constant 
hopping t we calculated analytically the frequency of the 
collective modes. It is of the form lu = \J (2t) 2 + lUq in 
the neutral system. In the case of intralayer pairing ui 2 , 
is proportional to the ratio — J/|Vji| while for interlayer 
pairing it is proportional to (Vj| — V± — 2J)/\V\\ + V±\ (in 
the limit that theses ratios are small) . The former result 
is in agreement with the result for the collective mode 
obtained by Wu and GriffinlI3 for t = 0. The latter result 
is new. In order to find an undamped collective mode 
it is important that its frequency is below the particle 
hole threshold \J (2t) 2 + (2A) 2 . As for realistic systems 
uiq > the best chances to observe a collective mode 
is for the case of intralayer pairing and small Josephson 
coupling J. However, one has to keep in mind that the 
Coulomb interaction between the layers leads to a further 
shift to higher frequencies. This Coulomb interaction is 
proportional to the distance between the layers but can 
be reduced by a large phononic polarizability of the in- 
termediate layers. 

In principle one has to distinguish two different modes 
which in the literature sometimes are called phase and 
amplitude mode. For u!q < the two collective modes, 
can both have frequencies below the particle-hole thresh- 
old for the same coupling parameters in the neutral sys- 
tem. However, in the case of constant t and equal values 
of the two gaps (pure intralayer or interlayer pairing) 
only one of these modes couples to the electronic polar- 
ization and will show up in the optical spectra. This is 
different for mixed intralayer and interlayer pairing and 
finite hopping (this occurs in our model for k-dependent 
hopping tk)- Then both phase and amplitude oscilla- 
tions couple to the charge oscillation between the layers. 
Therefore in the c-axis optical conductivity two collective 
mode peaks can appear below the particle-hole threshold 
for low enough Coulomb interaction. 

For a more realistic parameter region (negative J and 
strong Coulomb interaction) only one collective mode, 
the plasmon corresponding to charge fluctuations be- 
tween the layers, exists. In this case most of the weight 
of the polarization spectrum is concentrated in the plas- 
mon peak; the particle-hole excitations and in particular 
the peak at the threshold of particle-hole excitations are 
strongly suppressed. Therefore it is difficult to determine 
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the threshold energy for breaking up Cooper pairs or the 
superconducting gap with optical c-axis experiments. It 
is remarkable that the plasmon peak is much narrower 
in the normal state than in the superconducting state. 
That is due to the stronger quasiparticle dispersion in 
the superconducting state, allowing a much wider range 
of interband transition energies. 

Up to now the c-axis plasmon is not detected in the 
layered high-T c superconductors YBCO or BSCCO. One 
explanation could be, that in these compounds the single- 
particle hopping t is zero (suppressed by correlation ef- 
fects, as suggested by P.W. Anderson) and furthermore 
interlayer pairing dominates. The fluctuation of inter- 
layer Cooper-pairs does not produce density fluctuations 
between the layers. The weight of the polarization func- 
tion would then only be determined by the hopping 
matrix-element and therefore vanishes. But it is also 
possible (not investigated here) that the plasmon is over- 
damped by impurity scatteringo or more likely by in- 
elastic scattering effects at high frequencies due to anti- 
ferromagnetic spin fluctuations. 
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APPENDIX A: MATRIX ELEMENTS 

For the matrix Kij we have to calculate the following 
integrals: 

= ~ ^2 7; Tr {a-Jcr^G 1 {k, iu> n + iu s )ajG 2 {k, iuj n )a 3 } 

U) n k 

(Al) 

where <7$ are Pauli matrices and the matrices <7j are de- 
fined by: 

<jj = Uj for j = 0, 1, 3, a 2 = i(J2, &2 = —i&2 
Defining the integrals 

H^iius) = -5^5^[(iw n +iw,)iw n ±fi fc f 2fc ±AiA 2 ]/JV 

" w„ k 



(X/X)(iu s ) = — y^y^[(io;„ + iuj s )&k T iun£ik]/N 

" Un k 



(Z/Z)(iuj s ) = - Y^Yy^2 ± fcfcAJ/W (A2) 



with 

N = [(iw„ + iuj s ) 2 - E 2 lk ][{iuj n f 
The matrix K can be written as 



K 



E. 



2k \ 



( H++ 


-Y 


-z 


X 


Y 


-H-+ 


X 


z 


Z 


X 


-H— 


Y 


V x 


-Z 


-Y 


H+ 



(A3) 



/ 



The different signs come from the commutation relations 
for the Pauli matrices. Integrals of this type-are com- 
mon in vertex equations Jar superconductorsEa and for 
antiferromagnetic systemsa. 

The frequency summations are all of the general form 



J(k,iuj s ) = -3^2 



F(iu n + iuj s ,iui n ) 



f3 ^ {{iu> n + iuj s ) 2 - E 2 k )((iu n ) 2 



where the functions F(z\, z 2 ) are given by 



E 2 2k ) 



(A4) 



(X/X) 
(Y/Y) 
(Z/Z) 



F(zi,z 2 )= zxz 2 ± fut&ft ± AiA 2 

F(zi,z 2 )= z^2k T z 2 i\k 

F( Z1 ,Z2)= Zi A 2 t z 2 A ± (A5) 

F{z 1 ,z 2 )=€ lk A 2 ±&kAi 



After performing the frequency summations with help of 
Poisson's summation formula we obtain: 



J(k,iui s ) 



1 



,-. r , ( tanh 
8EikE 2 k \ 

F(— Eik, E 2 k) 



(3E 



lk 



tanh 



2 k 



2 2 

F{Eik 7 —E 2 k) 



E 



ik 



1 



tanh 



E 2k 
PE lk 



8Eik-E 2 k 

F(Eik,E 2 k) 



iuj s - Eik 

t ^$ E2k 
tann 



E- 



2k 



(A6) 



E 



lk 



E 



2 k 



F{-Ei kl —E 2 k) 
iui s + Ei k - E; 



2 k 



One notices that the first term containing contributions 
from the creation or destruction of two quasi-particles in 
different bands remains finite also for T — > 0. The second 
term with transfer of quasi-particles from one band to the 
other vanishes in that limit. Using the expressions for the 
function F(z', z) the different matrix-elements are easily 
calculated. One observes that the integrals H ±± ,Z, Z 
are even functions of uj s while the integrals X,X,Y,Y 
are odd functions. 
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APPENDIX B: VERTEX FUNCTIONS AND 
WARD IDENTITIES 

We derive some useful relations between different inte- 
grals of the matrix K, which are based on the conserving 
approximation used for the calculation of self-energies 
and vertex-functions. Similar relations for a one-band 
system are calculated e.g. in RetJ. 

The vertex equation in the neutral system in ladder 
approximation has the general form (|24|). 



T jl (iu s ) = D l1 -L{Y ij } 



(Bl) 



with 



VD 03 G(k, iuj n + iu s )T(iu s )G(k, ius n )D 03 
+VD 13 G(k, ioj n + iu; s )T(ioj s )G(k, iuj n )D 13 
+JD 33 G{k, iuj n + iu a )T(iu B )G(k, iuj n )D 33 (B2) 

We consider only vertex functions T, which are off- 
diagonal in the band indices: 



r = 



7 
7 



(B3) 



Then L is also off-diagonal in the band indices and can 
be written as 



L 



A 
A 



(B4) 



Now let us define a k-dependent vertex function: 

Tl(iuj s ) = G~ l {k, iuo n + iuj s )D 13 - D 13 G" 1 {k, iu n ) 

(B5) 

then the right upper corner is: 

7 £ = ioj s a 3 + 24do - t(Ai + A 2 )(T 2 (B6) 

where i k = (£ 2fc - 6fc)/2. 

Inserting T a into the r.h.s. of (B2) some of the Green's 
function cancel, and the integrals are reduced to those 
which also occur in the self-energies (|l4|). In particular 
we obtain for the 2x2 matrix in the right upper corner 

V = EE W (G 2 (k,iuj n )a 3 -a 3 G 1 (k,iuj n + iuj s )) 

k LO n 

+ V (Gi(fc, iuj n )a 3 - a 3 G 2 (k,iu> n + iuj s )) 



W 2 - V 2 

x { W [a 3 (-WE 2 + VSi) + (WEi - VZ 2 )<J 3 ] 
+V [tr 3 (-WSi + VY, 2 ) + (WZ 2 - V"Ei)<r 3 ] } 



W + V W - V 

-(A 2 + A x )ia 2 + — - — w (x 2 -x 1 )a (B7) 



W + V 



w-v 



On the other hand we can calculate L{T} in an other 
way. With the quantities B and B defined in (^7|) A is 
given by: 



A a = W'B + VB 



(B8) 



Because of the fc-dependence of the vertex function 7JJ 
we introduce the fc-dependent 4x4 matrix functions K% 
and K k , which are related to K and K by 

K = Y,K k , K = Y,Kk 



The same relations we needed for deriving the vertex 
equation (|4(]) holds for the /c-dependent quantities K k , 
Kk, 7JJ, and 7^. Therefore we find: 



\ a = (W'+gV)Y,K kl l 



(B9) 



where g is the diagonal matrix g — diag{—\, 1, 1, 1). Thus 
we arrive at the first identity: 

(W + gV) K kl a k = (x 2 -x 1 ,0,A 2 + A 1 ,0) t (BIO) 



Note on the l.h.s. stands W not W, because the factor 
W + gV cancels. 

In a similar way we find a second relation starting from 

T b k (iw s ) = iu n + iuj s )D w - D 10 G-\k, iuj n ) 

(Bll) 

Inserting this into the integral L we obtain 
A& = EE W '^G 2 (k, iLo n )a 3 - Va 3 Gt(k, 

k U>n 

-W'a 3 Gi(k,iuj n + iuj s )a 3 + Va 3 G 2 {k, iuj n + iui s )a 3 



w + v 
w-v 



(Si-X3a) 



(B12) 



or 



(B13) 



In the case of vanishing Josephson coupling and k- 
independent hopping t the first relation is identical with 
the Ward identity 



ioj s r p + 2itr j = r? 



(B14) 



which follows from the continuity equation 2tj = P (com- 
pare the discussion to (]§)). T p is here the polarization 
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vertex and T J the current vertex. The second relation is 
not a Ward identity based on a conservation law. It is 
valid only within the ladder approximation for the vertex 
equation. More relations are obtained by replacing the 
matrix D 13 by D 11 or D 12 . These relations, however, are 
not of practical use. 

From the above relations we obtain the following iden- 
tities between the different integrals: 



2H*++ 
2Y l 
2Z % 
2X* 



H 



++ 



10J S 



Yiu) s 
Ziui s 
Xiu> q 



Z2A|| 
X2A|| 
H— 2A|| 
F2A|| 

y2Aj_ 

X2A ± 
Z2Aj_ 



Xioj s 
Ziu> s 
Yiuj s 
H+-iLu s 
2X l 
2Z l 
2Y l 
2H t+ - 



{X2-X!)/W ± 



2A||/Wj| 



2A±/W± 


{X 2 -X 1 )/W ± 

(B15) 

where W\\ = Vy+J = W + V , W± = V±+J = W-V and 
A|| = (A 2 + Ai)/2, A_l = (A 2 - Ai)/2. Here the quan- 
tities K\j are integrals similar to Kij with the integrand 

multiplied by t^ = tk + (xi — x 2 )/2 before the momen- 
tum integration. In the case of momentum independent 
hopping matrix- element tk we have Kfj = tKij. 

The identities can be used to eliminate the badly con- 
verging integrals H h and H : 



H-+ = 



1 



H— = 



1 

WZ ~ A 2 - A! 
1 1 



W\\ Ai + A ; 



(iuj s Y + 2Z l ) 
(iuj s Y + 2Z t ) 



(B16) 
(B17) 



Here the integrals on the r.h.s. do not need a cut-off. 

In the two special cases of pure intralayer and inter- 
layer pairi ng one can derive from the definitions of the 
functions (|A5|) two further relations: 
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H 



-+ 



(A! + A 2 )Y/iuu s if Ax = A 2 (B18) 
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FIG. 4. Imaginary parts of the polarization function in the 
neutral (a) and charged (b) system at T = in the case of 
dominant intralayer pairing for different Josephson couplings 
J. The parameter sets (22,23) are used. The results are 
normalized by iVo = 13A?o. The different lines are shifted by 
equal units. The rhombs, crosses, and dashed lines in the 
figures correspond to various approximation formulas for the 
collective mode frequencies discussed in the text. 
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FIG. 5. Imaginary parts of the bare polarization function 
(a), the polarization function in the neutral system (b), and 
the real part of the optical conductivity in the charged sys- 
tem (c) in the case of dominant intralayer pairing for different 
temperatures T and positive Josephson coupling No J = 0.05. 
The parameters (22,23) are used. The different lines are 
shifted by equal units. In (c) we choose as distance between 
the layers d = 3A. 
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FIG. 6. Real part of the optical conductivity cr(w) in the 
case of dominant intralayer pairing for different temperatures 
T and negative Josephson coupling No J = —0.06. The pa- 
rameters (22,23) are used, d = 3A. The different lines are 
shifted by equal units. 
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perature T = 0. The parameters (22,23) are used. 
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FIG. 8. Imaginary part of the polarization function in the 
neutral (full lines) and charged (dashed lines) system in the 
case of dominant interlayer pairing for different temperatures 
and Josephson coupling J = 0. The parameters (22,23) are 
used. 
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